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I.  INTRODUCTION 


This  report  summarizes  the  work  done  on  Contract  No.  N00123-76-C-1371 

for  the  period  1  August  1976  to  31  July  1979. 

» 

The  overall  program  is  divided  into  two  areas.  These  are  1)  basic 
scattering  code  development,  and  2)  reflector  antenna  code  development. 

The  following  sections  present  background  material  and  describe  the  re¬ 
search  accomplished  in  each  of  these  two  areas. 

II.  PROGRAM  SCOPE 

The  scope  of  the  work  under  Contract  No.  N00123-76-C-1371  is  to  de¬ 
velop  the  necessary  theory,  algorithms  and  computer  codes  for  simulating 
antennas  at  UHF  and  above  in  a  complex  ship  environment.  A  milestone  chart 
for  the  three  year  program  is  shown  in  Table  1. 
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TABLE  I 

REVISED  MILESTONE  CHART  (04/12/78) 


III.  BACKGROUND 


The  performance  of  shipboard  communications  and  radar  systems  is  de¬ 
pendent  on  the  abilities  of  their  antennas  to  operate  effectively  in  a 
complex  electromagnetic  environment.  System  performance  at  UHF  and  above 
can  be  degraded  by  pattern  distortion  effects  due  to  blockage  and  scat¬ 
tering  by  shipboard  structures  or  by  strong  coupling  between  antennas [1]. 
In  many  cases,  these  detrimental  effects  can  be  minimized  by  the  judicious 
location  of  the  antennas.  This  task  is  complicated  by  the  large  number 
of  systems  that  are  used  on  a  modern  ship  where  many  systems  must  compete 
for  a  limited  number  of  prime  antenna  locations. 

The  two  approaches  that  can  be  used  for  modeling  shipboard  antennas 
at  UHF  and  above  are  scale  model  measurements  and  numerical  modeling. 

In  many  ways  these  are  complimentary  techniques.  Scale  model  measurements 
have  the  disadvantage  of  the  cost  and  time  required  for  a  thorough  design 
that  will  address  all  anticipated  problems.  It  is  very  desirable  to  fore¬ 
see  as  many  problems  as  possible  so  that  costly  remedial  solutions  can 
be  avoided.  Numerical  modeling,  on  the  other  hand,  is  a  relatively  cost 
effective  and  fast  way  of  modeling  the  many  possible  situations  that  can 
arise  in  a  shipboard  environment.  Measurements  need  be  made  only  when 
necessary  for  confirmation  of  results. 

Two  user  oriented  computer  codes  have  recently  been  developed  in 
order  to  allow  the  analysis  of  far  field  pattern  distortion  effects  of 
the  ship's  structure  on  an  antenna's  performance  at  UHF  and  above:  the 
NEC  Basic  Scattering  Code  and  the  NEC  Reflector  Antenna  Code.  The  develop 
ment  of  these  codes  has  been  accomplished  by  the  use  of  the  Geometrical 
Theory  of  Diffraction  (GTD)  [2,3,41 .  The  GTD  has  been  successfully  applied 
to  many  types  of  complex  modeling  problems  in  the  past[5,6,71.  It  is  a 
high  frequency  technique  that  allows  a  complicated  structure  to  be  approxi 
mated  by  basic  shapes  representing  canonical  problems  in  the  GTD.  These 
canonical  solutions  include  the  diffraction  at  flat  and  curved  wedges  and 
the  scattering  from  convex  curved  surfaces. 


IV.  BASIC  SCATTERING  CODE 


A.  Features  of  the  Code 

The  NEC  Basic  Scattering  Code  can  be  used  to  simulate  many  practical 
scattering  problems  by  using  basic  shapes  to  model  the  major  scattering 
features  of  a  complex  environment,  such  as  those  encountered  on  a  ship. 

The  basic  shapes  included  in  the  code  allow  the  analysis  of  scattering 
structures  that  can  be  modeled  by  perfectly  conducting  convex  and  concave 
flat  plate  structures,  and  a  finite  elliptic  cylinder.  The  basic  scat¬ 
tering  shapes  can  be  arbitrarily  positioned  in  space  to  form  the  various 
ship  structures  in  the  immediate  environment  of  the  antenna  under  analysis. 
For  example,  the  flat  plates  can  be  fashioned  into  box-like  shapes  that 
represent  the  ship's  superstructure,  platforms,  yardarms,  etc.  The  finite 
elliptic  cylinder  can  be  used  to  represent  smoke  stacks,  masts,  etc.  An 
example  is  illustrated  in  Figure  1. 

The  scattering  structures  are  assumed  to  be  large  in  terms  of  a  wave¬ 
length.  The  general  rule  is  that  the  lower  frequency  limit  of  this  solu¬ 
tion  is  dictated  by  the  spacing  between  the  various  scattering  centers 
and  their  overall  size.  In  practice,  this  means  that  the  smallest  dimen¬ 
sions  should  be  on  the  order  of  a  wavelength.  This  can  often  be  relaxed 
to  approximately  a  quarter-wavelength.  The  upper  limit  on  the  frequency 
is  dictated  by  the  fine  structural  details  that  can  be  modeled  by  the  basic 
scattering  shapes  available.  The  antenna  is  assumed  to  be  in  the  near 
field  of  the  scattering  structures,  but  the  scattering  centers  are  assumed 
to  be  in  the  far  field  of  the  antennas.  The  observer  is  assumed  to  be 
in  the  far  field  of  the  structures,  thus  the  scattered  rays  are  all  parallel 
leaving  the  structure  in  the  direction  of  the  observer. 

The  basic  scattering  code  is  capable  of  analyzing  a  wide  range  of 
antennas  since  it  is  basically  a  Green's  function  type  solution,  that  is, 
the  solution  has  been  developed  for  infinitely  small  current  elements. 

If  the  current  distribution  of  an  antenna  is  known,  such  as  from  the  NEC 


Moment  Method  Code[8l  ,  the  field  from  the  various  current  segments  can  be 
summed  to  give  the  total  field.  If  the  scattering  centers  are  in  the  far 
field  of  the  antenna  as  a  whole,  the  antenna  can  be  represented  by  its 
pattern  factor.  This  can  save  a  large  amount  of  computer  time  by  elim¬ 
inating  multiple  passes  through  the  diffraction  routines  for  each  infin¬ 
itesimal  current  element. 

There  are  two  documents  describing  the  NEC-Basic  Scattering  Code. 

Part  I  is  a  User's  Manual[9]  that  contains  a  detailed  description  of  the 
input  parameters  which  are  based  on  a  command  word  system,  an  interpre¬ 
tation  of  the  output,  and  example  problems.  The  example  problems  are  com¬ 
posed  of  sample  input  data  with  the  resulting  far  field  patterns  compared 
against  known  results  to  confirm  the  validity  of  the  code.  Most  users 
of  the  code  will  find  that  the  User's  Manual  is  sufficient  to  learn  how 
to  effectively  operate  the  code. 

The  Code  Manual  [10]  is  Part  II.  It  describes  the  FORTRAN  coding 
in  detail.  First,  however,  background  on  practical  aspects  of  the  GTD 
are  given  and  several  examples  are  shown  to  illustrate  how  the  various 
GTD  fields  superimpose  to  give  a  total  solution.  Next,  a  particular  GTD 
term  is  discussed  in  more  detail  to  show  the  general  concepts  involved 
throughout  the  code.  Then  the  code  manual  gives  an  overview  on  how  the 
code  is  organized.  It  describes  the  various  coordinate  systems  involved, 
how  a  general  subroutine  is  organized,  and  how  the  various  subroutines 
are  interrelated.  The  code  description  contains  for  each  subroutine: 

(1)  a  statement  of  purpose,  (2)  an  illustration  showing  the  geometry  in¬ 
volved,  (3)  a  brief  narrative  on  the  method  used,  (4)  a  flow  diagram,  (5) 
a  dictionary  of  major  variables,  (6)  a  listing  of  the  code.  Finally  the 
common  blocks  are  defined  and  system  library  functions  used  by  the  code 
are  listed. 

The  information  in  the  Code  Manual  will  be  of  primary  interest  to 
someone  attempting  to  modify  the  code.  It  will  also  be  helpful  when  the 
code  is  being  implemented  on  a  computer  system  on  which  the  coding  may 
not  be  compatible. 


The  basic  scattering  code  has  been  very  much  in  demand  during  its 
development.  The  agencies  to  which  various  generations  of  the  code  have 
been  sent  by  The  Ohio  State  University  are  listed  below.  The  first  gene¬ 
ration  of  the  code  was  sent  to  the  following: 

Westinghouse  (Baltimore,  Maryland) 

Electronic  Communications,  Inc.  (St.  Petersburg,  Florida) 

Naval  Research  Laboratory  (Washington,  D.C.). 

The  second  generation  of  the  code  was  sent  to  the  following: 

Martin  Marietta  Corp.  (Denver,  Colorado) 

Atlantic  Research  Corp.  (Rome,  New  York) 

H.  R.  B.  Singer  (State  College,  Pennsylvania) 

The  third  generation  of  the  code  was  sent  to  the  following: 

Naval  Post  Graduate  School  (Monterey,  California) 

Harris  Corp.  (Melbourne,  Florida) 

Boeing  Corp.  (Seattle,  Washington) 

General  Dynamics  (San  Diego,  California) 

Andrew  Corp.  (Chicago,  Illinois) 

Ohio  University  (Athens,  Ohio) 

Naval  Weapons  Center  (China  Lake,  California). 

B.  Example  Results 

The  NEC  Basic  Scattering  Code,  with  all  of  the  capabilities  mentioned 
above,  is  applicable  to  a  wide  range  of  practical  problems.  Of  course, 
with  such  a  general  code  not  all  possible  situations  can  be  anticipated 
a  priori.  Every  effort  has  been  made  to  include  the  dominant  scattering 
mechanisms  that  could  be  cost-effectively  added  to  the  code  to  give  a  good 
engineering  answer.  In  order  to  show  the  versatility  of  the  code  and  to 
test  its  accuracy  a  large  number  of  example  problems  were  run  and  compared 
with  measurements  over  a  wide  range  of  geometries.  Many  of  these  results 
are  included  in  the  User's  Manuals  [9].  The  code's  results  have  been  tested 
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against  measurements  made  at  OSU  and  whenever  possible,  they  have  been 
tested  against  measurements  made  at  other  facilities  or  against  results 
found  in  the  literature. 

For  example,  the  basic  scattering  code  has  been  checked  against  meas¬ 
urements  made  by  Bach  [11]  on  a  variety  of  satellite  shapes.  The  code's 
calculated  results  for  a  dipole  in  the  presence  of  an  eight  sided  box  as 
illustrated  in  Figure  2,  is  compared  against  measurements  in  Figure  3. 

The  calculated  results  for  various  locations  about  a  finite  circular  cylinder 
as  shown  in  Figures  4  a,  b,  and  c,  are  compared  against  measurements  in 
Figures  5  a,  b,  and  c,  respectively.  In  all  cases  the  agreement  is  very 
good,  confirming  the  validity  of  the  basic  scattering  code's  results.  The 
small  discrepencies  that  exist  between  the  patterns  can  be  contributed 
to  measurement  inaccuracies  and  to  the  lack  of  double  diffraction  effects 
in  both  the  geometries  and,  in  the  case  of  the  cylinder,  to  the  pseudo 
caustic  caused  by  the  coalescence  of  rays  on  the  endcap  rims. 


ik  * 


Figure  2.  Electric  dipole  in  the  presence  of  an  eight  sided  box. 
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electric  dipole  on  the  y-axis  parallel  to  the  x-axis  in  the  presence  of 
a  finite  circular  cylinder  (pattern  taken  in  y-z  plane). 


c  dipole  in  the  x-z  plane  parallel  to  the  x-axis  in  the  presence  of 
a  finite  circular  cylinder  (pattern  taken  in  the  y-z  plane). 


The  two  geometries  considered  above  are  composed  of  either  plates 
alone  or  a  finite  cylinder  by  itself.  In  order  to  test  a  problem  containing 
both  a  plate  and  a  finite  cylinder  a  series  of  measurements  were  made  at 
OSU  on  a  simplified  geometry  that  resembles  a  situation  that  may  arise 
on  shipboard.  The  geometry  is  shown  in  Figure  6.  This  is  a  very  difficult 
problem  because  of  the  multiple  bounces  of  the  fields  that  can  occur  be¬ 
tween  the  plate  and  the  cylinder.  This,  however,  illustrates  the  engineer¬ 
ing  solution  that  can  be  obtained  using  the  scattering  code.  For  comparison 
purposes,  the  separate  pieces  are  first  studied  individually.  The  cal¬ 
culated  and  measured  patterns  of  the  half-wavelength  dipole  are  compared 
in  Figure  7a.  The  results  for  the  square  plate  are  compared  against  measure¬ 
ments  in  Figure  7b.  The  results  for  the  finite  circular  cylinder  are  shown 
in  Figure  7c.  In  all  of  these  three  cases  the  comparisons  are  quite  good. 

The  radiation  pattern  of  the  dipole  in  the  presence  of  both  the  plate  and 
the  cylinder  is  compared  to  the  calculated  result  in  Figure  7d.  This  re¬ 
sult  contains  five  plate-cylinder  interaction  fields.  It  can  be  noted 
that  the  levels  are  good  and  that  the  fine  details  have  the  right  trends. 

This  verifies  that  the  code  can  provide  good  engineering  answers  to  rather 
difficult  problems. 

The  examples  shown  above  verify  the  ability  of  the  code  to  model  gen¬ 
eral  type  shapes.  The  following  examples  illustrate  how  the  code  can  be 
used  to  model  antennas  of  particular  interest  to  the  shipboard  antenna 
designer.  The  first  is  for  an  AS-2410/WSC-1  antenna.  The  geometry  used 
to  model  this  antenna  is  illustrated  in  Figure  8.  The  currents  on  the 
dipole  elements  of  the  antenna  were  obtained  from  the  NEC  moment  method 
code  by  Mr.  Jim  Logan  of  NOSC. 

First,  in  order  to  test  the  compatibility  of  the  two  codes  for  ex¬ 
actly  the  same  situation,  a  simpler  problem  was  used.  This  is  composed 
of  four  quarter  wavelength  dipoles  over  an  infinite  ground  plane.  The 


four  dipoles  are  shown  in  Figure  8  as  the  solid  lines  parallel  to  the  x- 
axis.  The  infinite  ground  plane  for  this  case  is  coincident  with  the 
square  plate.  The  dipoles  were  broken  into  a  total  of  20  segments  and 
the  frequency  was  0.2998  GHz.  The  far  zone  fields  obtained  from  the  moment 
method  code  and  the  fields  obtained  using  the  moment  method  current  weights 
in  the  basic  scattering  code  are  compared  in  Figures  9  a  and  b  for  the 
planes  <{>=0o  and  <f>=90°,  respectively.  The  basic  scattering  code  calculations 
for  the  directive  gain  is  within  0.2  dB  of  the  moment  method  results  con¬ 
firming  the  compatibility  of  the  two  codes. 

Next,  to  show  how  the  directive  gain  changes  when  the  four  dipoles 
are  placed  over  a  square  plate  over  an  infinite  ground  plane  as  shown  in 
Figure  8,  the  same  moment  method  current  weights  are  used  in  the  basic 
scattering  code.  It  is  assumed  that  the  current  distribution  obtained 
for  the  infinite  ground  plane  case  will  not  be  perturbed  appreciably  by 
slightly  changing  the  environment.  The  comparisons  between  the  infinite 
ground  plane  and  the  square  plate  cases  are  shown  in  Figures  10  a  and  b 
for  the  two  pattern  planes. 

For  the  case  of  the  crossed  dipoles  of  the  WSC-1  antenna  of  Figure 
8,  the  dipoles  were  broken  into  a  total  of  48  segments  at  a  frequency  of 
0.25  GHz.  The  far  zone  fields  for  the  crossed  dipoles  over  an  infinite 
ground  plane  obtained  from  the  moment  method  (AMP)  code  and  the  basic  scat¬ 
tering  code  are  compared  in  Figures  11  a  and  b.  Figure  11a  shows  the  ver¬ 
tical  component  and  Figure  lib  the  horizontal  component  of  the  circularly 
polarized  field.  Again  the  two  results  compared  to  within  0.2  dB.  The 
square  plate  results  are  also  shown  in  Figures  12  a  and  b  for  the  <j>=0° 
and  4>=45°  planes,  respectively.  Again  these  results  confirm  the  compati¬ 
bility  of  the  NEC  moment  method  and  basic  scattering  codes. 

A  few  measurements  for  the  WSC-1  antennas  are  available  from  NOSC 
[1].  The  simple  model  of  the  WSC-1  antenna  in  Figure  8  is  used  along  with 
an  analytic  representation  of  the  8  dipoles  to  compare  with  a  measured 


Figure  9a.  Comparison  of  the  directive  gain  of  four  dipoles  over  an 
infinite  ground  plane  obtained  from  the  NEC  moment  method 
code  and  from  the  basic  scattering  code  using  the  moment 
method  currents  (<j>=C  vertical  polarization). 
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Figure  10a.  Comparison  of  the  directive  gain  of  four  dipoles  over  an 
infinite  ground  with  four  dipoles  over  a  square  plate  over 
an  infinite  ground  plane  (<j>=0°,  vertical  polarization). 
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Figure  10b.  Comparison  of  the  directive  gain  of  four  dipoles  over  an 
infinite  ground  with  four  dipoles  over  a  square  plate  over  an 
infinite  ground  plane  (<t>=90°,  horizontal  polarization). 
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Figure  11a.  Comparison  of  the  directive  gain  of  four  sets  of  crossed 
dipoles  over  an  infinite  ground  and  over  a  square  plate 
over  an  infinite  ground  (<j>=0°,  vertical  polarization). 
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pattern  from  a  full  scale  model.  The  analytic  representation  of  the  dipole 
is  used  to  save  computer  time.  The  results  for  the  component  of  the 
circularly  polarized  field  are  compared  in  Figure  13.  The  frequency  is 
0.250  GHz  and  the  antenna  is  pointed  at  0=0°  (90°  elevation). 

This  same  model  is  used  to  compare  against  measured  results  from  a 
1/10  scale  model  [1].  The  results  for  the  circularly  polarized  field  (total 
field)  and  compared  in  Figures  14  a,  b  and  c  for  the  antenna  pointed  at 
0=0°,  30°,  and  60°  (90°,  60°  and  30°  elevation),  respectively.  The  full 
scale  frequency  is  0.320  GHz.  The  basic  feature  of  the  patterns  in  all 
of  these  comparisons  are  good.  However,  the  secondary  lobes  are  off  slightly. 
This  may  be  due  to  the  simplified  model  used  to  calculate  the  results  using 
the  basic  scattering  code  along  with  normal  experimental  errors.  For  example, 
it  is  difficult  to  simulate  circularly  polarized  antennas  because  the  low 
level  fields  depend  on  the  precise  phasing.  The  model  used  for  the  compu¬ 
tations  can  be  greatly  improved  if  more  details  are  known  about  the  experi¬ 
mental  models  and  set  up.  For  example,  the  polarization  of  the  pattern 
could  be  more  closely  matched  if  more  details  of  the  receiving  antenna 
are  known.  Also,  the  model  of  the  antenna  could  be  improved  by  adding 
more  detail,  such  as  using  a  circular  cylinder  to  represent  the  antenna 
pedestal  and  using  more  plates  to  model  the  back  plate,  etc.  These  capa¬ 
bilities  are  in  the  code  but  the  necessary  information  was  not  readily 
available. 

The  second  example  is  for  an  AS-1735/SRC  antenna.  The  geometry  used 
to  model  this  antenna  is  illustrated  in  Figure  15.  It  is  composed  of  four 
dipoles  mounted  around  a  mast  and  is  used  to  provide  an  omnidirectional 
pattern  in  azimuth  [12].  This  antenna  is  to  operate  over  the  frequency 
range  of  0.225-0.46  GHz.  This  means  that  for  a  mast  of  5.5  inches  in  radius, 
the  electrical  size  varies  from  0.10-0.19  wavelengths  in  radius  and  the 
dipoles  are  mounted  0.17-0.29  wavelengths  off  the  surface.  In  general. 
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Figure  13.  Comparison  of  the  measured  and  calculated  Eo  patterns  for  the 
ia/SC-1  antenna.  The  calculated  result  is  from  the  NEC-Basic 
Scattering  Code.  The  measurements  were  made  on  a  full 
scale  model  at  NOSC [1 ] .  The  antenna  is  pointed 
at  e=0  (90  elevation). 
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Figure  14a.  Comparison  of  the  measured  and  calculated  circularly 
polarized  pattern  for  the  WSC-1  antenna.  The  calculated 
result  is  from  the  NEC-Basic  Scattering  Code.  The 
measurements  were  made  on  a  1/10  scale  model  at 
NOSC [1 3-  The  antenna  is  pointed  at  0=0° 

(90  elevation). 
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Figure  14b.  Comparison  of  the  measured  and  calculated  circularly 
polarized  pattern  for  the  WSC-1  antenna.  The  calculated 
result  is  from  the  NEC-Basic  Scattering  Code.  The 
measurements  were  made  on  a  1/10  scale  model  at 
NOSCtlL  The  gntenna  pointed  at  e=30° 

(60°  elevation). 
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Figure  14c.  Comparison  of  the  measured  and  calculated  circularly 
polarized  pattern  for  the  WSC-1  antenna.  The  calculated 
result  is  from  the  NEC-Basic  Scattering  Code.  The 
measurements  were  made  on  a  1/10  scale  model  at 
N0SC[1].  The  antenna  is  pointed  at  0=60° 

(30  elevation). 


the  smallest  applicable  dimensions  for  an  asymptotic  solution  are  on  the 
order  of  one  wavelength.  However,  the  uniform  asymptotic  solutions  are 
often  used  for  sizes  on  the  order  of  a  quarter  of  a  wavelength  for  engi¬ 
neering  purposes.  In  order  to  test  the  basic  scattering  code's  usefulness 
for  this  problem,  the  far  zone  field  of  a  single  dipole  was  calculated 
in  the  azimuth  plane  at  the  lower,  middle  and  upper  ends  of  the  frequency 
range.  The  resulting  patterns  are  compared  against  the  exact  modal  solu¬ 
tion  in  Figures  16a,  b,  and  c  for  the  frequencies  of  0.225  GHz,  0.3125 
GHz,  and  0.400  GHz,  respectively.  These  results  produce  an  omnidirectional 
pattern  in  azimuth  for  the  four  dipole  case  as  is  expected.  Therefore, 
the  patterns  have  not  been  shown.  The  agreement  is  very  good  and  even 
though  the  asymptotic  results  should  be  viewed  with  caution  for  such  small 
dimensions,  these  results  seem  to  verify  the  usefulness  of  the  basic  scat¬ 
tering  code  in  giving  an  engineering  answer  to  this  problem. 

The  third  example  is  for  an  AS-390/SRC  antenna  mounted  on  the  yard 
arm  of  a  mast  as  shown  in  Figure  1.  Measurements  were  made  for  this  con¬ 
figuration  by  L.  S.  Hansen  of  NOSC  [1]  using  a  1/10  scale  model  at  a  fre¬ 
quency  of  4.0  GHz.  The  geometry  used  to  model  this  structure  for  the  basic 
scattering  code  is  illustrated  in  Figure  17.  The  AS-390  is  simply  modeled 
as  a  dipole  of  scale  model  length  1.925  inches  for  the  code.  The  far  zone 
field  calculated  in  the  azimuth  plane  for  the  antenna  in  the  presence  of 
a  circular  cylinder  representing  the  mast  is  compared  with  the  measured 
result  in  Figure  18.  The  results  calculated  for  an  antenna  in  the  presence 
of  the  circular  cylinder  and  a  flat  plate  representing  the  yard  arm  is 
compared  against  the  measured  result  in  Figure  19.  Note  that  the  result 
for  the  cylinder  alone  gives  the  basic  shape  to  the  pattern.  The  flat 
plate,  however,  appears  to  give  the  result  better  agreement  in  the  null 
depths  and  in  the  shape  of  the  lobe  at  180°.  The  lack  of  a  small  feature 
at  about  160  degrees  needs  more  study.  This  may  be  due  to  something  present 
in  the  measurements  that  is  not  modeled  by  the  code  or  it  could  be  a  higher 
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Figure  16b.  Comparison  of  the  exact  modal  solution  and  the  uniform 
asymptotic  solution  of  the  E0  pattern  in  the  azimuth  plane 
for  one  dipole  of  an  AS-1735  antenna.  The  frequency 
is  0.3125  GHz. 
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Figure  16c.  Comparison  of  the  exact  modal  solution  and  the  uniform 
asymptotic  solution  of  the  E0  pattern  in  the  azimuth  plane 
for  one  dipole  of  an  AS-1735  antenna.  The  frequency 

is  0.400  GHz. 
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Figure  18.  Comparison  of  the  measured  and  calculated  azimuth 
pattern  for  the  AS-390  antenna  example.  The  calculated 
result  is  from  the  NEC-Basic  Scattering  Code  including 
the  mast  only.  The  measurements  were  made  on  a  1/1.0 
scale  model  including  yardarms  at  NOSC 
by  L.  S.  Hansen  [1  ]. 


42 


/ 


-  MEASURED 

-  CALCULATED 


order  inl.or.irt  ion  between  the  plate  and  the  cylinder  that  has  not  been 
included  in  the  code.  It  is  also  possible  that  this  feature  is  due  to 
the  currents  produced  on  the  antenna  due  to  the  reflected  field  from  the 
cylinder  incident  on  the  antenna  structure.  This  could  be  studied  by  using 
the  NEC  moment  method  code  to  analyze  the  antenna  currents  in  more  detail 
and  then  to  use  those  currents  in  the  basic  scattering  code.  In  any  event 
the  agreement  in  Figure  19  is  good  for  most  engineering  applications. 

V.  REFLECTOR  ANTENNA  CODE  > 

A.  Features  of  the  Code 


The  reflector  antenna  code  has  the  capability  for  both  near  field 
and  far  field  computations  for  most  reflector  antennas.  One  important  feature 
of  the  code  is  its  capability  for  a  general  reflector  rim  shape.  Another 
important  feature  is  its  capability  to  input  a  practically  arbitrary  volu¬ 
metric  feed  pattern. 

The  computer  program  package  includes  a  user's  manual  and  a  code  manual. 
The  user's  manual  [13!  describes  the  operation  of  the  code  which  is  based 
on  a  command  word  system  to  input  data  into  the  code.  By  this  method  ad¬ 
ditional  computer  runs  can  be  made  by  making  small  changes  in  the  input 
data.  The  code  manual  [14]  describes  the  theory  used  in  the  code  and  docu¬ 
ments  a  detailed  explanation  of  the  code.  It  covers  the  main  program  and 
its  various  subroutines.  The  purpose  and  method  was  given  for  each  section 
of  the  code,  accompanied  by  a  flow  diagram,  a  key  variable  list  and  a  code 
listing  for  that  section. 

The  reflector  code  can  be  used  for  the  following  basic  applications: 

1.  Pattern  prediction  of  existing  reflector  antennas. 

?.  Reflector  antenna  design. 
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3.  Radiation  hazard  calculations. 

4.  EMC  or  coupling  calculations  with  small  antennas. 

The  far  field  capability  of  the  code  is  used  for  applications  1  and 
2  listed  above.  For  pattern  prediction  it  is  necessary  to  have  sufficient 
information  about  the  reflector  dimensions  and  feed  pattern. 

For  antenna  design  the  code  can  be  used  in  an  iterative  manner  to 
seek  a  practical  design  having  a  given  pattern  performance  goal  or,  the 
code  can  be  used  to  give  a  more  accurate  prediction  of  the  performance 
of  a  design  obtained  from  more  approximate  techniques. 

The  near  field  capability  of  the  code  can  be  used  to  EMC  and  radiation 
hazard  applications.  The  code  can  accurately  calculate  the  field  at  virtually 
any  point  that  is  at  least  one  diameter  from  the  reflector.  Since  the 
code  is  efficient  even  for  near  field  computations  it  eliminates  the  necessity 
to  rely  entirely  on  approximate  techniques  as  has  usually  been  done  in 
the  past. 

For  radiation  hazard  applications  the  code  is  used  to  calculate  the 
level  of  the  electric  field  or  the  power  density  at  the  near  field  point. 

For  EMC  or  coupling  calculations  the  power  density  incident  on  a  small 
antenna  is  first  calculated  using  the  code;  then  the  coupling  is  calculated 
by  multiplying  by  the  effective  aperture  of  the  small  antenna. 

B.  Computed  Far  Field  Pattern  Results 

The  results  of  the  general  reflector  code  for  several  example  antenna 
geometries  are  presented  in  this  section.  These  examples  illustrate  the 
versatility  of  the  code  for  computing  far  field  patterns  of  reflector  an¬ 
tennas  with  various  rim  shapes  and  feed  patterns. 
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The  first  example  is  the  common  case  of  a  circular  reflector  with 
an  on-axis  feed  as  shown  in  Figure  20.  A  GTD  analysis  was  previously  de¬ 
veloped  for  this  case  as  reported  in  Reference  15.  Calculated  results 
are  given  in  Reference  15  for  a  24"  diameter  reflector  with  an  F/D  ratio 
of  1/3  and  a  frequency  of  11  GHz.  The  far  field  patterns  as  computed  by 
the  general  reflector  code  are  shown  in  Figures  21  a  and  b  for  the  H-  and 
E-planes,  respectively.  The  results  from  the  general  reflector  code  were 
found  to  be  in  good  agreement  with  the  calculated  results  of  Reference 
15  without  aperture  blockage  as  shown  in  Figures  22  and  23.  Aperture  blockage 
and  feed  strut  scattering  effects  have  not  yet  been  included  in  the  general 
reflector  code.  The  comparison  of  the  computed  results  with  the  measured 
patterns  in  Figures  22  and  23  is  one  indication  of  the  validity  of  the 
general  reflector  code.  Other  comparisons  are  given  below  to  further  demon¬ 
strate  the  accuracy  of  the  general  reflector  code. 

The  second  example  is  an  offset  reflector  as  shown  in  Figure  24.  For¬ 
tunately,  sufficient  data  were  published  by  Silver  ri6i  to  completely 
specify  the  antenna  and  most  of  its  feed  pattern.  The  feed  pattern  data 
for  this  reflector  had  to  be  extracted  from  the  published  dB  contours  [161 
of  the  aperture  field.  The  feed  is  horizontally  polarized.  The  feed  patterns 
were  specified  by  using  the  option  for  an  offset  feed  tilt  angle  (PSIT=20°) 
and  the  feed  pattern  symmetry  option  ( I SYM=3 ) .  This  feed  symmetry  option 
specifies  the  feed  pattern  to  be  symmetrical  about  the  y-z  plane.  The 
E-  and  H-plane  patterns  as  computed  for  9.2  GHz  from  the  general  reflector 
code  are  shown  in  Figures  25  a  and  b.  The  corresponding  measured  patterns 
as  published  in  Silver  [161  are  shown  in  Figure  26.  The  measured  and  com¬ 
puted  patterns  are  in  good  agreement  considering  that  feed  strut  scattering 
effects  are  not  included  in  the  reflector  code.  The  full  180°  patterns 
as  computed  by  the  general  reflector  code  are  shown  in  Figures  27  a  and 
b.  The  feed  spillover  is  evident  from  the  lobes  near  90°  and  110°  in  the 
E-  and  H-plane,  respectively. 


46 


PHI  (DEG) 


Figure  21.  Far  field  patterns  of  circular  reflector  example 
computed  by  general  reflector  code. 
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Figure  22.  H-plane  pattern  of  a  parabolic  reflector  with  a  flanged 
waveguide  feed.  Computed  in  Reference  15. 


Figure  25.  Far  field  patterns  of  offset  reflector  example  computed  by 
general  reflector  code.  Main  beam  and  near  sidelobes. 
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Figure  26.  Measured  E-  and  H-plane  patterns  of  the  antenna 
shown  in  Figure  24  (from  Ref.  16). 


patterns  of  offset  reflector  example  computed  by  general  reflector 


The  next  example  is  the  offset  reflector  shown  in  Figure  28.  The 
projected  rim  shape  is  circular  as  shown,  although  the  rim  is  actually 
elliptical.  The  feed  is  horizontally  polarized  with  an  offset  angle  of 
51.5°.  A  specialized  computer  code  was  previously  developed  by  Mentzer 
et  al  [17]  which  is  applicable  only  for  this  particular  antenna.  Although 
the  mid-C  band  frequency  was  not  precisely  known,  several  frequencies  were 
tried  until  one  was  obtained  which  matched  the  beamwidths  of  the  patterns 
calculated  in  Reference  17.  The  patterns  computed  by  the  general  reflector 
code  are  shown  in  Figures  29  a  to  c.  They  are  in  excellent  agreement  with 
those  calculated  in  Reference  17  and  shown  in  Figures  30  a  to  c.  An  antenna 
gain  of  65  dB  was  computed  by  the  general  reflector  code;  and  the  gain 
of  the  feed  horn  was  computed  to  be  about  37  dB.  Thus  the  general  reflector 
code  predicts  the  antenna  gain  as  26  dB  above  the  feed  horn  gain.  This 
is  within  about  1  dB  of  that  shown  in  Figure  30. 
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Figure  29.  (Continued) 


The  depolarization  properties  of  offset  reflector  antennas  were  studied 
hy  Chu  and  Turrin  P81.  They  calculated  the  maximum  cross  polarization 
produced  by  a  linearly  polarized  feed  and  checked  their  result  with  a  meas¬ 
ured  cross  polarized  pattern  of  the  offset  reflector  with  an  offset  angle 
of  0c=45°  shown  in  Figure  31.  Kaufman  et  al  [19]  also  checked  this  case 
with  calculated  patterns  from  their  computer  program  for  circular  reflector 
antennas.  Their  computer  program  is  based  on  an  aperture  integration  analysis 
involving  the  quantizing  and  sorting  of  numerous  aperture  field  data. 

In  the  calculations  of  [18]  and  [19]  for  this  offset  reflector  example, 
a  circularly  symmetric  feed  pattern  was  assumed.  The  patterns  computed 
by  the  general  reflector  code,  as  shown  in  Figures  32  through  34,  are  in 
excellent  agreement  with  those  computed  in  [19]  and  with  the  measured  pat¬ 
terns  of  [181.  Note  that  excellent  agreement  is  achieved  for  the  cross 
polarized  beam  as  shown  in  Figure  33  for  the  4>=0°  plane.  There  is  no  cross 
oolarized  beam  in  the  <|>=90°  plane  for  a  symmetrical  feed  pattern. 


(a)  FRONT  VIEW  (b)  SIDE  VIEW 

Figure  31.  Offset  reflector  geometry  for  example  of  Chu  and 
Turrin  [18].  D  =  12  inches,  vIy=P£;=450,  f/D=0.25 
and  frequency  =  18. d  6Hz. 
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Figure  32.  Far  field  pattern  for  example  of  Reference  18, 
computed  by  general  reflector  code.  Principal 
polarization  is  <fr=0°  plane. 
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Figure  33.  Far  field  pattern  for  example  of  Reference  18, 
computed  by  general  reflector  code.  Cross  polarization 

in  4f0°  plane. 


Figure  34.  Far  field  pattern  for  example  of  Reference  18,  computed 
by  general  reflector  code.  Principal  polarization 
in  4>=90°  plane. 


No  cross  polarization  occurs  when  the  offset  reflector  is  fed  by  a 
circularly  polarized  feed;  but  the  depolarization  property  manifests  itself 
by  a  beam  shift  for  circular  polarization.  The  left  and  right  circularly 
polarized  beams  shift  in  opposite  directions  [18].  The  corresponding  cir¬ 
cularly  polarized  cases  as  computed  by  the  general  reflector  code  are  shown 
in  Figures  35  and  36.  The  calculated  beam  shift  in  each  case  was  0.35°, 
or  a  total  shift  of  0.70°  between  left  and  right  circular  polarizations. 
This  compares  well  with  the  value  of  0.75°  as  determined  from  measured 
patterns  [18].  The  latter  are  reproduced  in  Figure  37. 

C.  Computed  Near  Field  Results 

In  this  section  some  computed  near  field  results  from  the  general 
reflector  code  are  given.  In  some  regions  either  AI  or  GTD  can  be  used 
and  the  results  in  this  section  show  both  AI  and  GTD  calculations.  In 
addition,  some  near  field  results  obtained  by  the  Plane  Wave  Spectrum  (PWS) 
method[20]  are  available  for  the  case  of  a  uniform  circular  aperture. 

The  PWS  results  are  compared  with  those  of  a  circular  parabola  calculated 
by  the  reflector  code  using  a  very  large  F/D  ratio. 

The  large  F/D  example  consists  of  a  circular  reflector  with  D  =  10X 
and  an  isotropic  feed  located  at  a  focal  distance  F=10,000X.  The  near 
field  results  were  computed  for  constant  plane  cuts  at  Z  =  5.3,  10,  12, 

40X  as  shown  in  Figures  38  through  43.  Part  a  of  each  figure  shows  the 
near  field  as  calculated  by  AI  and  part  b  shows  that  for  GTD.  The  PWS 
results  for  the  uniform  circular  aperture  are  shown  in  Figure  44  for  the 
corresponding  Z  cuts.  As  seen  from  these  figures,  all  AI,  GTD  and  PWS 
are  in  excellent  agreement.  Note  that  the  minor  Ez  components  in  the  Z=10X 
(Figure  40)  and  Z=40X  (Figure  43)  cuts  also  check  with  that  of  PWS  (Figures 
44  b,  c).  Even  though  a  reflector  antenna  with  such  a  large  focal  distance 
is  not  practical  to  build,  this  example  does  show  the  validity  of  the  code 
for  near  field  calculations. 
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Figure  35.  Far  field  pattern  for  example  of  Reference  18, 
computed  by  general  reflector  code.  Left  circular 
polarization  in  <|>=0o  plane. 
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Figure  36.  Far  field  pattern  for  example  of  Reference  18 
computed  by  general  reflector  code.  Right  circular 
polarization  in  <J>=0°  plane. 
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Figure  38.  Near  field  component  E  for  10X  diameter 
circular  reflector.  Z^5.3X. 


72 


o 

cu 


I 


PHI  (DEG) 
D  (WL) 

GRID  (WL) 
Z  (WL) 

fll  ' 


0.0 

10.00 

0.500 

10.0 


UJ 


(a)  Aperture  Integration  (AI) 


RHO  (WAVELENGTHS) 

(b)  GTD 

Figure  39.  Near  field  component  E  for  10X  diameter 
circular  reflector.  2M0X. 
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Figure  41.  Near  field  component  E  for  1( 
circular  reflector.  ?=12x. 
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Figure  42.  Near  field  component  E  for  10x  diameter 
circular  reflector.  Z^40x. 
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Figure  43.  Near  field  component  E  for  10X  diameter 
circular  reflector.  /?=40A. 
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Figure  44.  PWS  calculations  for  near  field  of  uniform 
circular  aperture.  10X  diameter. 


The  next  example  presents  the  near  fields  of  the  24"  diameter  (22.35a) 
reflector  antenna  for  which  far  field  comparisons  are  given  in  Figures 
21  to  23.  This  reflector  has  an  F/D  ratio  of  1/3  and  the  patterns  correspond 
to  a  frequency  of  11  GHz.  Figure  45  shows  the  near  fields  of  this  reflector 
at  range  R=40A  by  AI  (a)  and  GTD  (b).  As  can  be  seen  the  AI  and  GTD  cal¬ 
culations  are  in  very  good  agreement.  The  reflector  code  would  normally 
use  GTD  over  the  whole  pattern  region  for  this  example  which  shows  the 
good  agreement  between  AI  and  GTD  even  in  the  near  axis  region.  Figure 
46  shows  the  near  fields  of  a  constant  range  at  R  =  100A.  The  agreement 
between  AI  and  GTD  is  also  good  in  the  AI  region  except  near  the  axis. 

Both  AI  and  GTD  would  normally  be  used  for  this  example:  AI  for  6<12° 
and  GTD  for  0>12°.  As  can  be  seen  from  Figure  46  AI  and  GTD  overlap  in 
a  wide  region. 
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(a)  Aperture  Integration  (AI) 


Figure  46.  Principal  near  field  component  for  22.35A  diameter 
circular  reflector.  R=100A. 
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